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PREFACE 


This  report  represents  research  carried  out  from  1 September  1975 
to  30  September  1976  under  AFOSR  Grant  76-2947,  "Investigation  of  the 
Propagation  Mechanism  of  Pressure  Fluctuations  in  Turbulent  Flow".  This 
contract  was  technically  monitored  by  Lt.  Colonel  R.C.  Smith  of  the 
Department  of  the  Air  Force,  AFOSR,  Bolling  Air  Force  Base,  Washington,  D.C. 

Ihe  authors  wish  to  acknowledge  the  cooperation  of  Grumman  Aerospace 
Corporation  for  allowing  us  the  use  of  their  Ambilog  200  facility  for 
data  conversion.  We  also  wish  to  thank  Mr.  T.  Keller  for  his  guidance 
and  assistance  in  preparing  the  necessary  programs. 


ABSTRACT 


A wind  tunnel  facility  (M  = 0.6  to  0.8)  with  extremely  low  noise 


New  York  University  Aerospace  and  Energetics  Laboratory.  The  facility 
is  to  be  used  for  the  purpose  of  investigating  the  propagation  mechanism 
of  pressure  fluctuations  in  turbulent  boundary  layers.  The  facility  was 
designed  as  an  induction  tunnel  in  order  to  minimize  the  influence  of 
the  downstream  region  of  the  wind  tunnel  on  the  flow  field  in  the  test 
section.  This  report  presents  a description  of  the  facility  and  the 


flow  quality  obtained,  and  a discussion  of  the  instrumentation  and 


perimental  procedures  being  used.  The  results  obtained  from  the 


analysis  of  several  measurements  are  also  presented.  This  includes 


spectra  of  the  velocity,  wall-pressure  and  wall-shear  fluctuations  and 


the  results  of  a digital  processing  scheme  which  conditionally  samples 


the  fluctuations  to  obtain  a description  of  turbulent  "bursts"  in  the 


boundary  layer 
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I.  INTRODUCTION 

Turbulent  flows  have  been  investigated  for  several  decades. 
Technically  speaking,  almost  every  important  fluid  flow  is  turbulent, 
and  such  flows  have  considerably  higher  surface  drag  than  laminar  flows. 
Experience  has  shown  that  turbulent  flow  noise  is  radiated  both  directly 
from  the  boundary  layer  and  indirectly  from  the  boundary  walls  which 
are  stimulated  by  the  turbulent  wall  pressure  fluctuations.  There  is  a 
relationship  between  turbulence  production,  the  occurrence  of  turbulent 
wall  pressure  fluctuations,  and  the  so-called  turbulent  shear  stress 
(Reynolds -stress) . By  influencing  the  turbulent  flows,  it  should  be 
possible  to  reduce  drag  and  turbulent  flow  noise.  Better  knowledge  of 
the  mechanism  of  turbulence  production  could  show  how  it  is  possible  to 
influence  a fully-developed  turbulent  flow,  for  example,  by  means  of 
suitably  designed  flexible  walls,  without  using  additional  energy. 

Before  any  research  effort  could  be  directed  towards  an  understand- 
ing of  the  mechanisms  by  which  drag  and  noise  reduction  could  be 
accomplished,  it  would  be  necessary  to  examine  carefully  the  structure 
of  the  boundary  layer  turbulence  under  normal  flow  conditions.  It  is 
the  intention  of  the  research  reported  here  to  primarily  establish  the 
details  of  the  structure  of  such  a boundary  layer. 

It  has  been  found  in  many  experiments  (Refs.  1-10)  that  large  scale 
disturbances  that  have  substantial  three-dimensional  characteristics 
called  "bursts",  appear  in  boundary  layer  flow.  The  "burst"  is 
described  as  a thick  low-speed  streak  lifting  up  from  the  surface. 

Such  a profile  is  unstable,  and  the  formation  is  probably  connected  with 
some  kind  of  instability  triggered  by  small  disturbances.  A basic 
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objective  of  this  investigation  is  to  study  the  mechanism  that  triggers 


such  a "burst"  and  what  effect  it  has  on  the  structure  of  the  boundary 


layer.  The  necessity  of  producing  real  two-dimensional  boundary  con 


ditions  is  an  important  point  of  the  experiment.  Many  so  called  two 


dimensional  investigations  are  made  on  flat  plates,  or  flat  tunnel 


walls.  In  these  experiments  the  presence  of  ends  or  corners  can  have 


substantial  effects  in  determining  some  of  the  properties  of  turbulence 


It  is  for  this  reason  that  an  axially  symmetric  configuration  has  been 


chosen  for  this  investigation 


Furthermore  it  has  been  determined  that  the  frequency  and  amplitude 


of  the  instabilities  associated  with  the  turbulent  "bursts"  can  be 


characterized  as  being  uniformly  distributed  through  the  flow.  These 


instabilities  have  been  shown  to  move  downstream  with  a constant  speed 


U , (i.e.,  the  flow  structure  repeats  itself  in  space  and  time).  No 


interaction  occurs  between  neighboring  "bursts"  and  U is  greater  than 


the  local  flow  velocity 


Black  (Ref.  10)  has  tried  to  make  an  energy  balance  by  choosing  a 


coordinate  system  based  on  one  of  the  burst  cells.  By  basing  his  energy 


balance  on  the  incoming  and  outgoing  primary  velocity  profiles  relative 


to  this  system,  he  performed  a mass,  momentum,  energy,  and  vorticity 


balance  to  deduce  the  vortex  flow  pattern  within  the  burst.  Based  on 


this  observation,  the  vortex  pattern  responsible  for  the  eruption  of 


fluid  from  the  sublayer  to  the  outside  flow  could  be  calculated,  and 


its  contribution  to  the  additional  shear  stress  within  the  boundary 


layer  found 


To  the  best  of  our  knowledge  the  research  which  has  been  performed 


up  to  now  in  investigating  these  "bursts 


has  involved  measurements 


made  either  at  the  wall  (Ref.  6)  or  in  the  boundary  layer  (Refs.  16-18). 

In  this  investigation  simultaneous  measurements  will  be  made  of  the 

pressure  and  shear  at  the  wall  and  of  the  velocity  in  the  boundary  layer. 

In  this  manner  it  will  be  feasible  to  get  a relationship  giving  the 

fluctuating  shear  stress  at  the  wall  in  terms  of  the  induced  local  pressure 

fluctuations.  The  tests  were  performed  in  a one  foot  tunnel  at  M = . 64 . 

00 

Present  day  experimental  research  in  turbulent  flows  relies  heavily 
on  electronic  measuring  devices,  and  therefore  any  report  which  pertains 
to  this  type  of  research  is  not  complete  unless  a discussion  of  the  in- 
strumentation being  used  is  included.  It  is  for  this  reason  that  a 
significant  portion  of  this  report  is  devoted  to  a detailed  discussion  of 
the  instrumentation  and  measurement  techniques  utilized  in  this  investiga- 
tion. 


II.  DESCRIPTION  OF  FACILITY 
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A.  WIND  TUNNEL 

1 . Wind  Tunnel  Set-up 

An  induction  wind  tunnel  was  constructed  for  the  present  experiment. 

It  was  built  oi  heavy  gauge  steel  material  The  inlet  section  of  the 
tunnel  consists  of  a circular  cone  5 ft  in  diameter  and  has  several  meshed 
screens.  The  cone  is  followed  by  a cylindrical  section  3 ft  in  diameter  in 
which  plastic  tubes  (1.5  inches  in  diameter)  arc  arranged  parallel  to  the 
center  line  of  the  wind  tunnel.  These  plastic  tubes  act  as  flow 
straighteners . Several  fine  me^h  screens  were  placed  downstream  of  the 
flow  straighteners  in  order  to  further  reduce  the  turbulence  scale.  A conical 
section  3 ft  in  length  connects  the  inlet  to  a pipe  with  an  inside  diameter 
of  12  inches  and  a length  of  30  ft  This  pipe  is  followed  by  a cylindrical 
test  section  which  is  5 ft  long,  has  an  inside  diameter  of  12  inches,  and  is 
fitted  with  four  rectangular  windows.  The  inner  surfaces  of  the  windows  were 
contoured  in  order  to  match  the  test  section  inside  diameter.  The  instru- 
mentation for  the  measurements  were  mounted  on  these  windows.  A three  degree 
half  angle  conical  diffuser  was  placed  downstream  of  the  test  section  to 

reduce  the  flow  velocity  prior  to  the  entrance  to  the  vacuum  sphere 
3 

(40,000  ft  ) The  entire  assembly  was  isolated  from  the  supporting  structure 
and  its  vibrations  were  damped  by  putting  thick  rubber  pads  at  the  supports. 
Sketches  of  the  wind  tunnel  and  test  section  are  given  in  Figs.  1 and  2 

In  order  to  minimize  the  noise  level  of  the  wind  tunnel,  the  flow 
passing  through  it  was  generated  by  induction  by  connecting  the  downstream 
end  of  the  tunnel  to  the  vacuum  sphere.  A 6°  half  angle  cone  was  put 
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inside  the  tunnel  immediately  downstream  of  the  test  section  to  create  a 
sonic  throat.  In  this  manner,  the  propagation  of  noise  from  sources  down- 
stream of  the  test  section  was  eliminated  at  the  throat. 
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To  start  the  wind  tunnel,  the  pressure  in  the  vacuum  tank  was  reduced 
to  a very  low  pressure.  A hydraulic  valve  located  between  the  vacuum  tank 
and  the  diffuser  was  opened  and  atmospheric  air  was  induced  into  the  wind 
tunnel.  Steady  flow  was  established  within  half  a second  after  the 
tunnel  was  started.  By  moving  the  conical  center  body  upstream  and  down- 
stream with  respect  to  the  test  section,  the  cross  section  of  the  sonic 
throat  can  be  changed  in  order  to  vary  the  free  stream  velocity  within  the 
test  section.  For  the  present  experiments,  a Mach  0.64  free  stream  was 
established  with  a vacuum  pressure  of  8 psia  and  the  smallest  contraction 
ratio  between  the  test  section  and  the  sonic  throat. 

2 . Wind  Tunnel  Calibration 

The  flow  in  the  test  section  vas  calibrated  with  a pitot  probe  and  a 
hot-film  probe  utilizing  a constant  temperature  anemometer.  The  free  stream 
Mach  number,  boundary  layer  profiles  of  the  mean  velocity  and  turbulence 
intensity  (u'/U)  in  the  test  region  were  measured.  Details  of  the  wind 
tunnel  calibration  are  presented  here. 

Pressure  orifices,  0.0312  inches  in  diameter,  have  been  flush  mounted 
in  the  streamwise  direction  along  the  inner  surface  of  the  test  section  to 
measure  the  surface  static  pressure.  These  measurements  were  used  in  order 
to  insure  that  no  appreciable  axial  pressure  gradient  exists  during  the 
test,  and  was  monitored  for  all  tests.  With  a wind  tunnel  total  pressure 
of  14.7  psia,  the  streamwise  surface  pressure  decreases  from  11.13  psia  to 
10.97  psia  in  a distance  of  3 ft.  The  free  stream  Mach  numbers  for  these 
tests  were  between  0.64  and  0.66. 
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In  order  to  determine  the  velocity  boundary  layer  thickness,  the  static 


and  total  pressure  distributions  across  the  boundary  layer  at  the  test 
section  were  measured.  A pitot  probe  was  traversed  in  a distance  of  6" 
normal  to  the  tunnel  inner  surface.  This  pitot  probe  was  connected  to 
pressure  transducers  calibrated  within  a range  of  0-15  psia  to  measure  the 
local  static  pressure  and  total  pressure.  Based  on  these  pressure  distri- 
butions, the  boundary  layer  Mach  number  profile  and  the  corresponding 
velocity  profile  with  the  assumption  of  constant  total  temperature  were 
calculated  and  shown  in  vig.  3. 

The  boundary  layer  turbulent  intensity  profile  has  been  measured  with  a 
cylindrical  alumina  coated  hot  film  probe  and  constant  temperature 
anemometer.  The  hot- film  probe  was  also  traversed  across  the  boundary 
layer.  The  DC  and  RMS  voltage  outputs  from  the  constant  temperature 
anemometer  were  also  measured  at  different  locations  normal  to  the  wall. 

The  local  turbulent  intensity  was  computed  according  to  King's  law.  The 
boundary  layer  profile  of  the  turbulent  intensity  is  presented  in  Fig.  4. 

A free  stream  turbulent  intensity  of  0.008  was  obtained  and  a maximum  tur- 
bulent intensity  of  0.08  has  been  found  near  the  wall.  The  free  stream 
turbulent  intensity  obtained  for  these  tests  were  to  the  authors'  knowledge, 
the  lowest  thus  far  obtained  at  these  velocities. 

3.  Flow  Quality 

Based  on  the  results  of  the  wind  tunnel  calibration,  the  velocity 
boundary  layer  thickness,  displacement  thickness,  and  momentum  thickness 
have  been  computed.  The  boundary  layer  velocity  profile  could  also  be 
described  by  a 1/7  power  law.  The  parameters  of  the  boundary  layer  flow 
in  the  test  section  are  presented  in  the  following  table: 


* 


u 

00 


0 


Re0 

uj/v 

u'/U 


B . INS  TRUMENTATION 


675  ft/sec 
0.334  ft 
0.047  ft 
0.032  ft 
1.06  x 105 
1.61  x 10^  1/sec 
0.008 


It  has  been  found  in  many  experimental  investigations  that  large  scale 
quasi-ordered  disturbances  (called  "bursts")  having  substantial  three- 
dimensional  characteristics  appear  in  boundary  layer  flow.  Because  of  the 
intermittent  and  coherent  nature  of  the  bursts,  usual  time  and  space 
averaged  measurements  tend  to  "smear  out"  their  essential  features.  It 
has  become  necessary  to  utilize  various  visual  (Refs.  12-15)  and  condi- 
tional sampling  (Refs.  16-19)  techniques  in  order  to  reveal  more  clearly 
the  "bursting"  characteristics  of  the  turbulent  boundary  layer.  However, 
the  previous  investigations  involved  measurements  made  either  at  the  wall 
or  in  the  boundary  layer.  To  our  knowledge,  this  will  be  the  first  time 
that  simultaneous  measurements  of  the  fluctuating  pressure  and  shear  at 
the  wall,  and  fluctuating  velocity  in  the  boundary  layer  are  made  in 
order  to  study  the  generation  and  structure  of  turbulent  "bursts". 

1.  Data  Acquisition 

The  general  set-up  of  the  data  recording  instrumentation  is  shown  in 
Fig.  5.  Basically,  two  systems  are  required  - one  for  measuring  and 
recording  the  velocity  in  the  boundary  layer  and  the  fluctuating  shear  on 
the  wall  of  the  tunnel  and  another  for  measuring  and  recording  the 
fluctuating  pressure  on  the  wall. 
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a.  Velocity  and  Shear 

The  mean  and  fluctuating  components  of  the  streamwise  velocity  in  the 
flow  are  measured  with  hot-film  probes  and  constant  temperature  anemometers. 
Two  types  of  probes  are  available,  those  with  an  alumina  coated  cylindrical 
sensor  and  one  with  an  alumina  coated  wedge  type  sensor.  In  the  former 
case,  miniature  probes  (Thermo -Systems  type  TSI-1260)  are  available  for  use 
on  a rake  holding  five  probes  (see  Figs.  6 and  7).  Two  rakes  have  been 
constructed  - one  spanning  a distance  of  2 inches  from  the  wall  of  the 
tunnel  and  the  other  1 inch  as  shown  in  Fig.  6.  The  constant  temperature 
anemometer  systems  (DISA  and  Thermo-Systems)  used  in  conjunction  with  these 
probes  all  include  a linearizer  and  signal  conditioner  (amplifier/filter). 
The  DC  components  of  the  linearized  signals  from  the  hot-film  probes  are 
recorded  on  a 30  channel  Honeywell  1612  visicorder  and  are  monitored  during 
each  test  on  several  meters.  The  visicorder  is  also  used  to  record  the 
static  pressure  along  the  wall  of  the  test  section  to  insure  that  the  pro- 
per running  condition  is  achieved  in  each  test.  The  fluctuating  component 
of  each  of  the  signals  is  amplified  and  filtered,  and  then  recorded  on  one 
of  7 channels  of  a Phillips  Ana -Log  7 magnetic  tape  recorder.  The  RMS 
levels  of  the  fluctuating  signals  are  also  monitored  during  each  test  on 
in-line  meters. 

The  fluctuating  shear  on  the  tunnel  wall  is  measured  with  a flush 
mounted  hot-film  sensor  (Thermo -Systems  type  TSI-1237)  as  shown  in  Fig.  6. 
This  measurement  is  made  by  replacing  one  of  the  hot-film  probes  in  Fig.  5 
with  this  flush  hot-film  sensor  since  the  instrumentation  required  is  the 
same  as  that  described  in  the  previous  paragraph  for  the  hot-film  probes. 
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b.  Pressure 


Flush-mounted  piezoresis tive  pressure  transducers  (Kulite  type  MIC- 
080-5)  are  used  to  measure  the  pressure  fluctuations  on  the  tunnel  wall. 

The  diameter  of  the  pressure  sensitive  area  of  the  transducers  is  .040" 
and  they  are  rated  for  a pressure  differential  of  5 psi.  A natural 
frequency  of  200  kHz  is  quoted  by  the  manufacturer.  Five  such  transducers 
are  available  and  can  be  placed  at  various  spacings  in  a cross-shaped 
pattern  on  a plug  in  the  tunnel  wall  below  the  hot-film  probes  (see  Figs.  6 
and  7).  (This  plug  is  also  used  to  mount  the  flush  hot-film  sensor).  The 
signal  from  each  of  the  transducers  is  put  through  an  amplifier/filter  and 
then  recorded  on  one  of  7 channels  of  a Honeywell  5600  tape  recorder. 

In-line  RMS  meters  allow  the  signal  levels  to  be  monitored  during  each  test. 
The  input  power  for  the  transducers  is  provided  by  batteries  which  can 
supply  up  to  9 volts  DC. 

The  two  magnetic  tape  recorders  have  the  capability  of  recording  both 
in  FM  or  Direct  mode  at  speeds  as  high  as  60  ips . All  the  elements  of  the 
instrumentation  have  been  chosen  so  as  to  allow  sensing  and  recording  of 
fluctuations  at  fr  mencies  as  high  as  40,000  Hz  (ojS  /U  ==  18).  The 

CO 

synchronization  of  the  data  on  the  two  recorders  is  accomplished  by 
simultaneously  recording  two  signals  on  each  recorder,  a 100  kHz  sine 
and  a pulse  train  with  a period  of  from  1 to  3 sec.  The  sine  wave  is  pro- 
vided by  a Hewlett-Packard  3310A  function  generator  and  the  pulse  train  by 
a General  Radio  1340  pulse  generator  (see  Fig.  5.) 

2 . Data  Analysis 

The  tape  recorded  data  can  be  analyzed  with  the  systems  shown  sche- 
matically in  Fig.  8.  Spectral  analyses  of  the  pressure,  shear  and  velocity 
fluctuations  are  obtained  using  a General  Radio  model  1900  wave  analyzer 
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and  model  1521  graphic  level  recorder,  while  futo  and  cross-correlations  are 
performed  with  a Saicor  SAI-42  correlation  and  probability  analyzer.  A 
continuous  display  of  5 simultaneous  data  channels  is  provided  by  an  ink 
plotter  which  by  slowing  down  the  tape  recorder  sufficiently  on  playback 
can  yield  a resolution  of  up  to  3 msec/inch.  The  sixth  channel  of  the 
plotter  is  used  to  display  the  recorded  pulses  which  serve  as  reference  time 
marks . 

The  most  important  aspect  of  the  data  analysis  involves  the  digitization 
of  the  recorded  fluctuations  and  the  subsequent  manipulation  of  this  data 
on  a large  digital  computer.  To  accomplish  this  an  Ambilog  200  computer  with 
an  analog-to-digital  converter  is  utilized  to  obtain  discrete  time  sequences 
of  the  fluctuations.  The  digitization  process  is  triggered  by  one  of  the 
recorded  pulses  (chosen  depending  on  the  portion  of  test  to  be  analyzed); 
and  the  rate  at  which  samples  are  taken  can  be  established  in  one  of  two 
ways,  either  by  a timing  loop  within  the  digitizing  program  itself  or  by 
using  the  recorded  100  kHz  sine  as  a clock.  When  it  is  necessary  to 
preserve  the  simultaneity  of  the  data  on  two  channels  digitized  at  different 
times  or  of  the  data  on  the  two  tape  recorders,  the  latter  procedure  is 
chosen.  (The  need  for  these  procedures  will  be  discussed  in  greater  detail 
in  the  next  section).  With  a sufficiently  large  ratio  of  record-to-playback 
speed,  it  is  possible  to  achieve  sampling  rates  as  high  as  200,000  samples/sec. 
Up  to  24,000  samples  can  be  made  and  stored  during  one  pass  through  the 
digitizing  program.  All  digital  data  is  transferred  to  magnetic  tape  which 
serves  as  the  input  medium  to  a CDC  6600  computer  where  final  numerical 
analysis  of  the  data  is  performed. 

The  overall  experimental  set-up  including  the  instrumented  tunnel  test 
section  and  data  acquisition  instrumentation,  as  well  as  some  of  the  data 
analysis  equipment  are  shown  in  Figs.  9 and  10. 


10 


III.  EXPERIMENTAL  PROCEDURES 


A.  CALIBRATION  OF  INSTRUMENTATION 

1.  Hot-Film  Probes  and  Constant  Temperature  Anemometers 

The  principle  of  the  hot-film  probe  and  anemometer  is  based  on  the 
convective  heat  loss  from  an  electrically  heated  film  (or  wire)  caused  by 
the  flow  of  gas  or  liquid  surrounding  the  wire.  Normally,  the  finite 
thermal  capacitance  of  the  transducer  will  affect  the  sensitivity  and  the 
phase  relation  so  that  changes  in  flow  velocity  can  be  measured  only  up  to 
a frequency  of  around  1000  Hz.  The  constant  temperature  anemometer  employs 
a technique  that  minimizes  the  effect  of  thermal  lag,  thereby  increasing  the 
upper  frequency  limit  by  a factor  of  approximately  100.  The  hot-film  probe 
uses  a conducting  film  on  a ceramic  substance^r  a quartz  rod  with  a 
platinum  film  on  the  surface.  The  hot-film  probe  sensor  has  a large  dia- 
meter and  is  particularly  helpful  in  reducing  the  breakage  of  the  sensor  due 
to  the  impact  of  the  flow.  For  the  present  experiment,  TSI-1260  miniature 
hot-film  probes  with  0.002  inch  diameter  by  0.04  inch  long  cylindrical 
sensors  were  used.  This  probe  was  capable  of  measuring  velocity  up  to 
1000  ft/sec  and  was  operated  with  an  overheat  ratio  of  1.5. 

A DISA  type  55D41  calibration  wind  tunnel  was  used  to  calibrate  the 
miniature  probes.  This  tunnel  gives  a maximum  velocity  of  700  ft/sec  with 
an  optimum  turbulent  intensity  of  0.003  at  300  ft/sec  and  100  kHz  bandwidth. 
The  hot-film  sensor  was  placed  perpendicular  to  the  mean  flow  direction  and 
the  mean  flow  velocity  was  measured  by  means  of  a pitot  tube. 

First  the  frequency  response  of  the  bridge  output  circuits  of  the 
anemometers  was  adjusted.  A square  wave  generator  and  oscilloscope  were 
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connected  to  the  anemometer  for  accurate  adjustment  so  that  no  damped 
oscillations  occur  on  the  output  signal  and  an  optimum  probe  time  constant 
was  obtained.  The  anemometers  upper  frequency  limit,  in  conjunction  with 
the  probe  being  used,  is  determined  by  this  time  constant.  An  upper 
frequency  limit  of  approximately  20  kHz  was  obtained  for  the  miniature 
probes  used  in  the  present  experiment. 

After  the  anemometers  frequency  response  has  been  adjusted,  the 
bridge  output  was  connected  to  a linearizer  to  obtain  a linear  relation 
between  the  mean  velocity  and  DC  voltage  output.  The  linearizer  output 
was  also  bypassed  to  an  RMS  meter  to  measure  the  root -mean-square  value 
of  the  output  signal  in  order  to  compute  the  turbulent  intensity.  These 
calibration  procedures  are  the  same  as  those  described  in  detail  in 
Ref.  20  where  a full  description  of  the  characteristics  of  the  flush  hot- 
film  sensor  is  also  given. 

2 . Pressure  Transducers 

The  Kulite  transducers  to  be  used  for  this  investigation  are  checked 
with  both  a static  and  dynamic  calibration.  In  the  static  calibration, 
the  pressure  differential  between  the  transducer  reference  tube  and  the 
sensing  face  was  varied  from  zero  to  the  rated  pressure  and  the  output 
voltage  noted  at  several  points.  The  output  voltage  was  found  to  vary 
linearly  with  the  applied  pressure  and  the  applied  voltage  to  the  trans- 
ducer bridge  circuit.  The  five  MIC-080-5  transducers  were  found  to  have 
a sensitivity  on  the  order  of  0.60  mV/psi/V(input) . One  MIC-080-25 
transducer  which  is  available  was  also  calibrated  and  its  sensitivity 
was  found  to  be  0.45  mV/psi/V. 
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The  dynamic  calibration  consisted  of  utilizing  a Bruel  & Kjaer 
pistonphone  which  supplied  an  oscillating  pressure  having  a frequency  of 
approximately  250  Hz  and  a precise  level  of  124  dB.  The  RMS  level  of  the 


output  voltage  from  the  transducer  exposed  to  this  pressure  field  was 
measured  and  yielded  a sensitivity  of  142  dB  below  lV/|ibar  for  the 
MIC-080-5  transducers  and  145  dB  below  lV/|jbar  for  the  MIC-080-25  trans- 
ducer. Both  figures  are  given  for  an  input  Voltage  of  9 volts  DC  and 
are  equivalent  to  the  static  sensitivities  listed  above.  The  frequency 
response  of  the  transducers  was  not  checked,  but  the  specified  natural 
frequency  of  200,000  Hz  would  insure  a flat  frequency  response  up  to 
40,000  Hz  (uj5*/U  = 18). 

B.  PROCEDURES  FOR  MEASUREMENT  OF  TURBULENT  FLUCTUATIONS 

After  the  various  instrumentation  systems  were  calibrated,  preliminary 
recordings  were  made  of  the  fluctuating  pressure  and  shear  on  the  tunnel 
wall  and  of  the  fluctuating  velocity  in  the  boundary  layer.  Recordings 
were  made  in  both  FM  and  Direct  mode  and  at  the  highest  speed  possible  on 
the  two  recorders  so  that  a data  bandwidth  of  at  least  DC-100,000  Hz  was 
available.  Except  for  the  DC  component  which  was  filtered  before  re- 
cording, the  amplifier  bandpasses  chosen  were  within  this  range.  It  was 
necessary  to  have  a wide  data  bandwidth  during  these  first  tests  so  that 
a good  indication  could  be  obtained  of  the  entire  spectral  range  of  the 
fluctuating  quantities,  in  particular  of  their  high  frequency  extent. 

Since  recordings  eventually  had  to  be  made  in  FM  mode  (to  allow  playback 
at  lower  speeds  without  distortion)  it  was  necessary  to  know  whether  the 
limited  high  frequency  response  in  the  FM  mode  (i.e.,  up  to  10  kHz  on 
the  Phillips  recorder  and  40  kHz  on  the  Honeywell)  would  be  sufficient  to 
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include  all  the  components  of  interest  in  the  fluctuations.  This  informa- 
tion was  also  valuable  in  determining  the  rate  at  which  the  data  would  have 
to  be  digitized. 


The  recordings  were  analyzed  to  obtain  the  overall  RMS  level  and 
spectral  distributions  of  the  fluctuations.  RMS  meter  readings  taken  both 
during  the  recordings  and  then  on  playback  were  used  in  conjunction  with 
calibration  data  to  obtain  the  overall  RMS  levels  of  the  various  fluctuat- 
ing quantities,  while  the  wave  analyzer  system  described  in  Section  II 
was  utilized  to  obtain  the  spectral  distribution  of  the  fluctuations.  A 
bandwidth  of  50  Hz  was  normally  used  for  the  spectral  analyses,  although  a 
10  Hz  bandwidth  was  also  employed  for  some  low  frequency  analyses.  The  FM 
recordings  were  analyzed  to  give  the  spectra  for  low  frequencies  while 
the  Direct  recordings  gave  corresponding  information  at  the  high  frequency 
end.  These  results  were  combined  to  yield  the  spectral  distributions  of 
the  fluctuating  quantities  over  a rather  broad  range. 

These  preliminary  recordings  were  also  displayed  on  a two  channel 
storage  oscilloscope  to  provide  qualitative  information  on  the  time  histories 
of  the  fluctuating  quantities,  and  to  verify  the  feasibility  of  syncroniz- 
ing  two  data  channels  on  different  recorders.  This  was  done  for  the  entire 
range  of  playback  speeds  to  insure  that  no  distortion  occurred  to  the 
recorded  signals  when  played  back  at  slower  speeds.  For  signals  recorded 
in  FM  mode  no  significant  distortion  was  found  even  for  ratios  of  record- 
to-playback  speeds  as  high  as  32  or  64.  As  was  expected  however,  the 
Direct  recordings  did  show  increasing  low  frequency  distortion  with  in- 
creasing ratios  of  record -to-playback  speed.  To  help  identify  the  im- 
portance of  various  frequency  ranges  on  selected  events  in  the  time 
histories  of  the  fluctuating  quantities,  the  recorded  fluctuations  were 
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also  passed  through  a bandpass  filter  before  being  displayed  on  the  oscil 


These  time  averaged  analyses  did  not  require  that  the  various  measure 


ments  be  made  simultaneously.  But  in  order  to  obtain  spatial  as  well  as 


temporal  resolution  of  the  "bursts"  in  turbulent  boundary  layers  it  is 


necessary  that  simultaneous  measurements  be  made  at  several  points  in  the 


flow.  The  following  capabilities  have  been  developed  for  simultaneous 


measurements  in  this  investigation:  up  to  five  wall  pressure  fluctuation 


measurements  together  with  either  five  streamwise  velocity  measurements  in 


the  boundary  layer  or  four  velocity  measurements  and  one  shear  fluctuation 


measurement  on  the  tunnel  wall 


Once  the  required  set  of  transducers  and  probes  are  placed  in  position 


in  the  tunnel  test  section  certain  procedures  are  followed  for  each  test 


First  the  signal  conditioning  elements  of  the  instrumentation  are  each 


checked  with  a test  signal.  The  filter  bandpasses  on  the  signal  conditioners 


are  normally  set  at  2-20,000  Hz  (u)6  /U 


0009-9.)  for  the  velocity  and  shear 


013-13.)  for  the  pressure  fluctua 


tions.  Just  before  flow  is  initiated  by  opening  a valve  to  the  vacuum 


sphere,  the  three  recorders  (the  two  magnetic  tape  recorders  and  the  visi 


corder)  are  started  and  the  anemometers  switched  to  operating  condition 


At  about  the  time  when  steady  flow  conditions  are  reached  the  pulse  generator 


is  turned  on  and  begins  to  mark  a channel  on  each  tape  recorder  with  a pulse 


every  three  seconds.  A 100  kHz  sine  is  also  being  simultaneously  recorded 


on  one  channel  of  each  recorder  at  all  times.  Up  to  ten  fluctuating 


quantities  (five  on  each  tape  recorder)  and  five  DC  anemometer  outputs  (on 


the  visicorder)  are  recorded  during  any  one  test.  During  the  test  itself 


meter  readings  are  also  taken  of  the  DC  anemometer  outputs  and  the  RMS  levels 


of  the  fluctuating  quantities.  An  entire  test  lasts  a total  of  approxi- 
mately 40  seconds  with  steady  flow  conditions  being  maintained  for  about  15 
seconds.  The  steady  portion  of  the  test  is  identified  with  the  use  of  one 
of  the  visicorder  traces  which  records  the  static  pressure  along  the  test 
section  wall  during  the  entire  test. 

To  date,  successful  recordings  have  been  made  for  the  following  test 
configurations : 

Four  hot-film  probes  on  the  1"  rake  and  a pressure  transducer 
flush  mounted  on  the  wall  directly  below  the  probes. 

Three  hot-film  probes  on  the  1"  rake,  a flush  hot-film 
sensor  directly  below  the  probes,  and  a pressure  trans- 
ducer approximately  .5"  upstream. 

- Five  pressure  transducers  in  a cross-shaped  pattern 
(equally  spaced  .15"  apart,  see  Fig.  6)  without  any 
rakes  in  the  flow. 

Five  hot-filra  probes  on  the  2"  rake,  and  the  five  pressure 
transducers  of  the  previous  configuration  centered  directly 
below  the  probes. 

Preparations  for  several  other  tests  are  in  progress  and  the  possibility 
of  repeating  some  tests  for  a lower  flow  velocity  (U  c&  100  ft/sec)  is 

CO 

being  considered. 

The  analysis  of  the. a tests  requires  the  digitization  of  the  recorded 
fluctuations  into  discrete  time  sequences.  This  is  accomplished  utilizing 
the  Ambilog  200  facility  described  in  Section  II.  The  five  data  channels 
and  the  pulse  channel  from  one  of  the  recorders  are  first  simultaneously 
displayed  on  a six  channel  ink-plotter.  This  is  done  for  both  tape  re- 
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corders  and  for  sections  of  recordings  falling  within  the  period  of  steady 
flow  conditions.  The  digitization  itself  required  that  a choice  be  made 
among  several  possible  procedures.  This  involved  choosing  the  way  in  which 
the  data  sampling  was  to  be  timed  and  whether  only  one  or  several  data 
channels  would  be  sampled  during  a single  execution  of  the  digitizing  pro- 
gram. One  of  the  first  procedures  to  be  tried  consisted  of  sampling  a 
single  data  channel  with  the  sampling  interval  being  timed  either  by  an 
external  clock  or  by  a timing  loop  within  the  program  itself.  Although 
initial  synchronization  between  successively  sampled  data  channels  was 
provided  by  having  the  same  recorded  pulse  serve  as  a trigger  for  the 
digitization  process,  it  became  apparent  that  because  of  the  lack  of 
synchronization  between  the  timing  mechanism  and  the  recorded  data  any 
real  time  unsteadiness  such  as  in  tape  speed  would  eventually  result  in 
the  loss  of  simultaneity  between  the  data  channels.  One  possible  solution 
to  this  problem  was  to  have  the  program  digitize  all  five  data  channels  on 
one  recorder  at  the  same  time.  But  since  a finite  time  interval  was  re- 
quired by  the  program  between  sequential  samples  of  each  channel  (result- 
ing in  a time  skew  between  the  data  from  the  five  channels)  and  because 
the  problem  described  above  would  still  exist  between  data  on  separate 
recorders,  this  procedure  was  not  implemented. 

It  was  finally  decided  that  the  most  acceptable  solution  would  be  to 
record  a 100  kHz  sine  on  each  recorder  to  serve  as  a clock  for  timing  the 
sampling  interval.  Of  course  the  key  to  this  procedure  is  that  the  sine 
wave,  like  the  pulses  which  serve  as  coarse  time  marks,  must  be  recorded 
on  both  tape  recorders  at  the  same  time  the  data  is.  In  this  way,  it 
provides  an  absolute  tape-time  reference  for  all  ten  data  channels  on  the 
two  recorders.  The  data  on  any  one  channel  can  be  digitized  using  this 


' 


17 


clock  in  the  following  way.  One  data  channel  and  the  pulse  and  100  kHz  sine 
channels  from  a given  recorder  are  fed  into  the  analog-to-digital  converter. 
The  program  in  core  begins  by  sampling  the  pulse  channel  and  comparing  the 
the  sampled  value  to  a pre-set  threshold.  When  the  threshold  is  exceeded 
the  pulse  has  been  reached  and  the  program  then  begins  to  sample  the  100  kHz 
sine.  At  every  sampled  zero  crossing,  or  at  any  multiple  number  of  zero 
crossings,  of  the  sine  wave  the  program  jumps  to  the  data  channel,  samples 
a value  and  stores  it  in  a core,  and  then  returns  to  sampling  the  sine 


wave.  The  program  continues  in  this  way  for  a pre-set  number  of  samples  of 
the  data  channel  up  to  a maximum  of  24,000.  To  obtain  another  set  of 
samples  of  the  same  channel,  the  procedure  is  repeated  with  a delay  in  the 
program  from  the  time  the  pulse  is  detected  to  when  sampling  of  the  data 
is  begun.  This  is  accomplished  by  counting  an  appropriate  number  of  zero 
crossings  of  the  sine  wave  before  the  first  sample  of  data  is  taken.  This 
entire  procedure  is  then  repeated  for  every  data  channel. 

The  pulse  and  sine  wave  channels  are  sampled  at  a rate  determined 
solely  by  the  time  required  to  perform  certain  program  operations.  Since 
this  time  is  much  larger  than  the  10  |asec  period  of  the  sine  wave,  it  is 
necessary  in  the  above  procedure  to  slow  down  the  tape  recorders  suffi- 
ciently on  playback  to  allow  the  required  sampling  resolution  of  the  sine 
wave.  For  this  reason,  ratios  of  record-to-playback  speeds  as  high  as  32 
or  64  must  be  used.  A maximum  data  sampling  rate  of  200,000  samples/sec 
would  be  possible  under  these  circumstances.  This  would  provide  five 
samples  within  one  cycle  of  40,000  Hz  (u)6  /U  = 18)  components  of  the  data. 

Once  the  digital  data  is  stored  on  magnetic  tape,  it  can  be  read  on 
a CDC  6600  computer,  converted  to  CDC  format,  and  numerically  analyzed 
with  Fortran  programs.  Basically,  a conditional  processing  scheme  similar 
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to  one  devised  by  Kaplan  and  Laufer  (Ref.  21)  is  utilized  to  infer  the 
occurrence  of  a "burst"  in  the  data.  This  involves  calculation  of  the 
variable  interval  time  average  (VITA)  variance  (see  Ref.  17)  defined  as: 


VAR  , (t) 


_ a'2(t)-ru'(t)r 


u 


,2 


where 


and 


t+T/2 

u’ (t)  = ^ f u'  (T)dT 
T Jt-T/2 


(T  = total  data  sample  time) 


From  this  a detector  function  is  calculated  according  to: 

1 for  VAR  . (t)  > V 
. u'  o 

IVARu,(t)  = { 

0 for  VAR  , (t)  < V 
u o 


(1) 
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where  V^  is  a preselected  threshold  level.  For  a physically  meaningfull 
average  time,  T,  and  threshold  level,  V , periods  of  non-zero  IVAR  should 
indicate  the  presence  of  "bursts"  in  the  fluctuating  quantity.  The  detector 
function  could  then  be  used  to  determine  the  mean  period  between  bursts, 
and  a correlation  of  the  variance  for  several  simultaneous  ^pressure  or 
streamwise  velocity  fluctuation  measurements  would  yield  an  estimate  of 
the  speed  of  propagation  of  bursts  along  the  boundary  layer  wall  or  normal 
to  the  wall.  An  analysis  of  instantaneous  velocity,  pressure  and  shear 
fluctuations  during  a burst  could  reveal  useful  information  on  the  mechanism 
of  burst  generation  and  on  the  structure  of  bursts. 
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IV.  RESULTS  AND  DISCUSSION 


A.  TIME  AVERAGED  ANALYSES 

The  recorded  data  from  the  various  tests  can  be  analyzed  in  two  basic 
ways.  Conventional  time  averaged,  analyses  (RMS  level,  spectral  density, 
auto  and  cross-correlation)  can  be  performed  for  comparison  with  similar 
work  done  by  other  investigators.  From  this  it  can  be  determined  if  the 
proper  instrumentation  systems  have  been  chosen  for  the  measurements,  if 
typical  turbulent  flow  conditions  have  been  established,  and  what  the  over- 
all boundary  layer  turbulence  characteristics  are.  To  obtain  any  informa- 
tion about  the  quasi-ordered  structure  of  the  turbulence  it  is  necessary  to 
analyze  the  data  in  a way  that  will  reveal  any  intermittency  or  coherence 
in  the  turbulent  fluctuations.  This  required  the  application  of  new  condi- 
tional sampling  schemes  to  the  data. 

Preliminary  recordings  of  the  fluctuating  pressure  and  shear  on  the 
tunnel  wall  and  of  the  fluctuating  velocity  in  the  boundary  layer  have  been 
analyzed  to  obtain  overall  RMS  levels  and  spectral  distributions.  Typical 
power  spectra  for  the  velocity,  shear  and  pressure  fluctuations  are  shown 
in  Figs.  11-14,  15  and  16,  respectively.  The  similarity  between  the  spectra 
for  the  shear  and  velocity  fluctuations  is  evident,  the  energy  falling-off 
after  approximately  uj6  /U  = .3  at  a rate  of  about  6 dB/octave.  Some  un- 

CO 

certainty  exists  in  the  accuracy  of  the  spectra  at  high  frequencies  due  to 
the  difficulty  in  optimizing  and  matching  the  frequency  response  charac- 
teristics of  the  various  anemometers.  This  is  particularly  evident  for  the 
spectrum  of  the  velocity  fluctuations  at  y/6  = .14  (Fig.  12)  which  shows 
that  the  anemometer  in  this  case  did  not  provide  as  high  a frequency  response 
as  the  others  A comparison  of  these  results  for  the  u-f luctuation  spectra 
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with  those  of  Serafini  (Ref.  11)  shows  good  overall  agreement. 

The  spectrum  of  Fig.  16  shows  that  the  pressure  fluctuations  contain 
significantly  more  high  frequency  energy  than  either  the  velocity  or  shear 
fluctuations.  Although  the  fall-off  in  energy  is  also  at  a rate  of  about 
6 dB/octave,  this  does  not  occur  until  around  uo6  /U  =3.  It  has  been  found 
by  several  investigators  (Refs.  5,  6,  and  11)  that  transducer  size  has  a 
significant  effect  on  the  resolution  of  wall  pressure  fluctuations.  A com- 
parison of  various  measurements  seems  to  indicate  tnat  if  the  transducer 
diameter  is  too  large,  a loss  of  resolution  of  small-scale  pressure  fluctua- 
tions will  occur.  This  manifests  itself  as  consistently  lower  measured 

•Jf 

spectral  densities  at  high  frequencies,  that  is,  above  ufi  /U  a»  1.  Emmerling 

(Ref.  6)  has  summarized  the  available  results  in  terms  of  the  nondimens ional 

parameter  du  /v,  where  d is  the  transducer  diameter.  For  values  of  this 
T 

parameter  above  approximately  150  the  overall  RMS  level  of  the  wall  pressure 

2 

fluctuations  is  always  measured  to  be  around  0.005q  (q  = '2  p U ).  Below 

CO  CO  CO  CO 

this  value,  the  measured  RMS  level  increases  linearly  as  du  /v  is  lowered  - 

levels  as  high  as  .Olq  have  been  measured  for  du  /V  “=  50.  This  increase  is 

oo  T 

attributed  to  the  increased  resolution  of  intense  small-scale  fluctuations, 
or  equivalently  of  high  frequency  spectral  components,  by  the  smaller 
transducers . 

For  the  present  investigation  the  condition  du^/v  *=  50  would  require  a 

transducer  diameter  on  the  order  of  .007".  Since  a transducer  this  small 

could  not  be  obtained,  it  was  decided  upon  consideration  of  sensitivity  and 

durability  requirements  to  use  Kulite  MIC-080-5  transducers  whose  pressure 

sensitive  area  has  a diameter  of  .040".  Although  this  gives  a value  of 

du  /v  equal  to  approximately  300,  it  was  planned  that  if  resolution  of  the 
T 

small-scale  fluctuations  was  required  to  achieve  the  goals  of  this 


investigation  a pinhole  surface  port  could  be  used  to  reduce  the  effective 
transducer  diameter.  The  main  obstacle  to  using  a surface  port  is  that 
the  pinhole-transducer  cavity  system  must  be  designed  so  that  its  resonant 
frequency  is  significantly  higher  than  any  frequency  of  interest  for  the 
pressure  fluctuations.  The  feasibility  of  designing  a system  which  meets 
this  requirement  for  our  flow  conditions  is  being  investigated.  The  results 
presented  in  Fig.  16  (as  well  as  in  the  rest  of  this  report)  are  for  the 
transducer  mounted  flush  with  the  tunnel  wall.  A comparison  of  this  wall 
pressure  spectrum  with  those  reviewed  by  Willmarth  (Ref.  5)  shows  that  some 
resolution  of  high  frequencies  is  being  lost  by  our  transducer.  The 
relatively  low  RMS  level  of  .00580^  also  bears  out  this  fact. 

The  spectral  analyses  showed  that  the  FM  mode  on  the  Phillips  tape 
recorder  was  adequate  for  recording  velocity  and  s!«^  “ fluctuations,  and 
that  the  FM  mode  on  the  Honeywell  tape  recorder  (with  its  wider  bandwidth) 
would  have  to  be  used  to  record  the  wall  pressure  fluctuations.  The  spectral 
content  of  the  data  also  indicated  that  the  digitization  would  have  to  be 
done  at  a sampling  rate  of  50,000  samples/sec  or  greater  in  order  to  obtain 
an  adequate  resolution  of  the  highest  frequencies  in  the  fluctuations. 

B.  CONDITIONAL  SAMPLING  ANALYSIS 


The  most  important  aspect  of  this  investigation  involves  the  conditional 
sampling  of  the  turbulent  fluctuations  to  obtain  a description  of  turbulent 
"bursts".  The  value  of  the  investigation  lies  not  only  in  the  fact  that 
simultaneous  measurements  are  made  of  the  velocity,  wall-pressure,  and  wall- 
shear  fluctuations  to  aid  in  obtaining  a better  description,  but  also  in 


the  range  of  flow  conditions  studied.  Table  I lists  two  typical  investiga- 
tions: one  by  Blackwelder  and  Kaplan  (Ref.  17)  of  conditionally  sampled 


? 
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I 

\ 
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turbulent  velocity  profiles  and  another  by  Emmerling  (Ref.  6)  of  the  in- 
stantaneous distribution  of  wall-pressure  fluctuations  under  a turbulent 
boundary  layer.  In  comparison,  the  present  investigation  is  of  a boundary 
layer  with  a much  higher  free  stream  velocity  (U  = 600-800  ft/sec)  and 

Reynolds  number  (R  s*  10"*),  thus  extending  the  range  over  which  the  "bursty" 
0 

character  of  turbulent  boundary  layers  has  been  studied. 

One  characteristic  of  interest  is  the  mean  period  between  "bursts". 

Lu  and  Willmarth  (Ref.  16)  have  found  that  a consistent  estimate  of  the 

- * 

average  time  interval  between  relatively  large  "bursts"  is  TU^/6  ^ 32 

(or,  equivalenty,  TU  /6  5 as  expressed  by  Laufer  and  Badri  Narayanan 

OO 

(Ref.  22)).  One  column  of  Table  I presents  the  ratio  of  this  period  to 
the  data  sampling  interval  and  shows  that  for  the  sampling  rates  possible 
in  the  present  study  an  accurate  temporal  resolution  of  the  "bursts" 
would  be  obtained. 

To  date,  the  data  from  one  of  the  test  configurations  listed  in 
Section  III-B  have  been  completely  digitized  and  analyzed.  This  was  the 
case  of  four  stre<-.nwise  velocity  measurements  within  a distance  of  1" 
from  the  tunnel  wall  and  one  pressure  fluctuation  measurement  on  the 
tunnel  wall  directly  below  the  velocity  measurements.  The  positions  of 
the  velocity  measurements  and  the  mean  velocity  and  overall  RMS  level 
of  fluctuations  measured  at  each  point  were  as  follows: 


y/6 

.025 

.08 

.14 

.25 

u/u 

CO 

.62 

.71 

.76 

.82 

y u'2 

.073 

.053 

.044 

.033 

U 


The  velocity  fluctuations  were  measured  within  the  frequency  range 


.009-4.5)  and  the  pressure  fluctuations  within  the 
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013-13).  The  overall  RMS  level  of  the 


pressure  fluctuations  was  approximately  / p'  / U 


A short  sample  of  the  traces  obtained  on  the  ink  plotter  directly 


from  the  tape  recorder  is  shown  in  Figs.  17a  and  b (the  two  figures 


actually  represent  six  simultaneous  traces).  The  duration  of  the  sample 
is  approximately  250  msec.  Since  the  pressure  fluctuations  were  recorded 
on  both  tape  recorders  for  this  test,  the  five  data  traces  shown  are  from 
one  recorder.  The  sixth  trace  is  of  the  pulse  channel  and  shows  one  of 


the  pulses  used  as  reference  time  marks.  Qualitatively  there  is  a great 


degree  of  correlation  in  the  time  histories  of  the  four  velocity  fluctua 


tion  measurements,  while  correlation  between  the  pressure  and  velocity 


Once  data  is  available  from  the  five  simulta 


neous  pressure  fluctuation  measurements  and  from  the  shear  fluctuation 


correlation  between  characteristic  events  in 


the  velocity  fluctuations  and  changes  occurring  at  the  wall  (e.g.,  in  the 
pressure  gradient  fluctuations)  will  be  more  readily  apparent. 

The  recorded  data  was  digitized  with  a sampling  interval  At  = 25  pisec 


(AtU  /fi  = .357  or  At  u /v  = 40).  Approximately  16,000  samples  (400  msec) 


were  taken  of  each  data  channel.  Two  sections  of  the  digitized  data  are 


plotted  in  Figs.  18a-h.  (Note:  the  first  four  figures  are  from  tU  /6  = 0 


to  1000  and  the  last  four  from  tU  /6  = 4700  to  5700).  The  plots  were 


made  on  the  CDC  computer  once  the  fluctuating  quantities  were  normalized 


with  respect  to  their  maximum  value.  For  the  velocity  fluctuations  at 


025,  .08,  and  .14,  the  vertical  scales  are  about  the  same  with  + 1 


representing  approximately  + 110  ft/sec.  For  y/6  = .25,  + 1 represents 


+ 75  ft/sec.  The  scale  for  the  pressure  fluctuations  is  such  that  + 1 is 


equivalent  to  + .075  psi 


The  results  obtained  from  the  digital  processing  analysis  are  shown  in 


4 


Fig.  19-23.  Each  figure  consists  of  eight  plots,  the  first  four  (a-d)  are 

* 

for  the  time  interval  5 < tU  /6  < 1005  and  the  last  four  (e-h)  for 

00 

4705  < tU  /6  < 5705.  The  results  have  been  presented  for  these  two  time 

00 

intervals  because  they  are  representative  of  those  obtained  for  the  entire 

test.  The  uppermost  trace  in  each  figure  is  of  the  fluctuating  quantity 

which  has  been  normalized  with  respect  to  its  maximum  value  over  the  entire 

16,000  samples.  Below  this  is  shown  the  variable  interval  time  average 

variance  (VAR)  of  the  fluctuating  quantity  as  defined  by  Eq.  1 of  Section 

III-B.  This  too  has  been  normalized  with  respect  to  its  maximum  value  so 

that  its  amplitude  varies  between  zero  and  + 1.  The  averaging  time,  T,  in 

Eq.  1 was  chosen  to  be  TU  /6  = 3.57  for  the  results  presented  here. 

00 

Several  other  values  were  tried  before  it  was  decided  that  this  averaging 
time  gave  the  most  consistent  results  for  the  variance.  The  lowermost 
trace  is  of  the  detector  function  (IVAR)  associated  with  the  VITA  variance 
through  the  relationship  expressed  in  Eq.  2 of  Section  III-B. 

It  can  be  seen  that  large  values  of  VAR(t)  indicate  steep  changes  in 
the  fluctuating  quantity,  that  is,  strong  streamwise  accelerations  in  the 
case  of  the  velocity  fluctuations  (Figs.  19-22)  and  high  temporal  pressure 
gradients  for  the  wall-pressure  fluctuations  (Fig.  23).  Since  various 
studies  have  reported  that  large  and  rapidly  changing  velocity  and  pressure 
fluctuations  are  associated  with  turbulent  "bursts",  it  is  commonly 
assumed  that  a period  of  large  amplitude  for  VAR(t)  indicates  the  presence 
of  a "burst"  in  the  fluctuating  quantity.  Once  a reasonable  decision  is 
made  as  to  what  constitutes  a large  amplitude  (that  is,  once  a threshold 
level,  Vq,  is  chosen),  a detector  function  can  be  defined  as  in  Eq.2. 

This  function  then  serves  to  mark  the  time  when  a "burst"  is  deemed  to 
occur.  For  the  results  presented  in  Figs.  19-23  the  threshold  was  chosen 
so  as  to  yield  the  best  correlation  between  the  detector  function  of  the 
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five  fluctuating  quantities. 

The  detector  function  can  be  used  to  obtain  the  mean  period  between 

"bursts"  for  each  of  the  fluctuating  quantities.  For  the  wall  pressure 

.*  * 

fluctuations  this  period  was  found  to  be  approximately  TU  / 6 = 32.5, 

CO 

— ft 

while  for  the  velocity  fluctuations  the  period  ranged  between  TU  /6  =34 

00 

— ft 

at  y/6  = .025  and  TU  / 6 = 36  at  y/6  = .25.  These  results  are  in  good 

CO 

agreement  with  those  reported  by  other  investigators.  For  example, 

_ * 

Emmerling  (Ref.  6)  has  found  a mean  period  of  TU  / 6 = 27  for  the 

00 

"occurrence  of  characteristic  wall  pressure  structures",  and  Lu  and 

- .'k 

Willmarth  (Ref,  16)  have  reported  a mean  period  of  TU  /o  = 32  for  rela- 

00 

tively  large  "bursts"  in  their  velocity  fluctuation  measurements. 

Although  it  is  evident  that  a significant  correlation  exists  between 
the  "bursting"  events  detected  in  the  four  velocity  fluctuation  measure- 
ments, a thorough  description  of  the  propagation  and  structure  of  the 
"bursts"  cannot  be  attempted  until  the  analysis  of  other  measurements 
included  in  the  full  test  schedule  (see  Section  III-B)  is  completed.  This 
will  include,  of  course,  a test  with  five  simultaneous  wall  pressure 
fluctuation  measurements  and  another  where  the  wall  shear  fluctuations 
are  also  measured.  These  measurements  will  be  invaluable  in  determining 
the  correlation  between  the  "bursts"  found  in  the  velocity  fluctuations 
along  the  boundary  layer  and  those  detected  on  the  wall.  The  simultaneous 
pressure  measurements  will  also  allow  the  gradient  of  the  wall  pressure 
fluctuations  to  be  included  in  the  analysis. 

In  addition  to  the  completion  of  the  full  test  schedule  and  digital 
analysis  of  the  collected  data,  the  following  investigations  are  in 
progress.  Several  of  the  tests  will  be  repeated  at  a significantly  lower 
free  stream  velocity  (U  r*  100ft/sec)  to  determine  how  any  observed 
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characteristics  of  the  "bursts"  scale  with  velocity.  This  will  also  serve 
to  check  any  inadequacies  (e.g.,  at  high  frequencies)  of  the  instrumenta- 
tion. Along  these  lines,  the  possibility  of  designing  a pinhole  surface 
port  to  decrease  the  effective  diameter  of  the  pressure  transducers  is 
being  investigated.  Since  it  is  critical  for  the  analyses  being  performed 
that  the  simultaneity  of  the  measured  fluctuations  be  preserved,  it  was 
necessary  to  determine  if  the  phase  characteristics  of  the  individual  in- 
strumentation systems  were  the  same.  All  the  elements,  except  for  the 
anemometers,  were  checked  and  found  to  have  no  relative  phase  shift.  The 
phase  characteristics  of  the  anemometers  are  in  the  process  of  being 
investigated  using  a procedure  described  in  Ref.  23. 

The  results  presented  in  this  report  are  preliminary  in  nature  but 
clearly  demonstrate  that  the  capability  has  been  established  for  examining 
the  structure  of  the  turbulent  boundary  layer  under  consideration.  The 
flow  has  been  made  truly  two-dimensional  and  the  free-stream  turbulence 
has  been  reduced  to  a relatively  low  level.  The  elaborate  instrumentation 
systems  required  have  been  set-up,  calibrated  and  shown  to  give  an 
acceptable  representation  of  the  turbulent  fluctuations.  Most  importantly, 
the  procedures  and  programs  necessary  to  perform  the  critical  conditional 
sampling  analyses  on  a digital  computer  have  been  established  and  ex- 
ercised. Once  the  full  test  schedule  is  completed  and  all  the  acquired 
data  analyzed  with  these  programs,  the  structure  of  the  turbulent  "bursts" 
can  be  better  described  and  their  generating  mechanism  more  clearly 
understood.  With  this  information  about  the  boundary  layer  turbulence 
characteristics  clearly  established,  it  will  then  be  possible  to  proceed 
with  the  next  step  in  the  investigation  which  will  be  to  attempt  to  in- 
fluence the  characteristics  by  imposing  external  disturbances. 


V.  CONCLUSIONS 


A wind  tunnel  facility  (M  = 0.6  to  0.8)  with  extremely  low  noise 


levels  (u'/u  = .08%)  has  been  completed  and  tested.  This  facility  is 


to  be  used  for  the  investigation  of  the  propagation  mechanism  of 


fluctuations  in  turbulent  flow.  Results  of  the  calibration  of 


the  facility  are  included  in  the  report  and  indicate  the  excellent 


applicability  of  the  facility  to  this  investigation.  A successful  in 


s trumentation  package  has  been  developed  and  tested  for  the  measurements 


of  the  velocity,  shear  and  pressure  fluctuations  in  the  turbulent  boundary 


layer.  Procedures  have  also  been  developed  to  produce  discrete  time 


sequences  of  the  fluctuations  which  are  conditionally  sampled  to  obtain 


the  frequency,  amplitude  and  structure  of  the  instabilities  associated 


with  turbulent  "bursts" 


Initial  tests  have  been  performed  at  M = 0.64  and  indicate  that  the 


spectra  of  the  turbulent  fluctuations  and  the  mean  period  between  "bursts 
is  in  agreement  with  results  obtained  by  other  investigations  at  lower 


free  stream  velocities.  These  previous  results  were  obtained  from  measure 


ments  made  either  at  the  wall  or  in  the  boundary  layer.  Our  conclusions 


are  based  on  simultaneous  measurements  of  the  pressure  at  the  wall  and 


the  velocity  in  the  boundary  layer.  Additional  tests  have  been  completed 
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8/ll(=k8*) 


4 


ltaneous  Traces  of  Velocity  and  Wall  Pressure  Fluctuations  Obtained  from  Analog  Recording 


Digitized  Velocity  and  Wall  Pressure  Fluctuations 


Digitized  Velocity  and  Wall  Pressure  Fluctuations  (500  < — — < 750) 


Pressure  Fluctuations  (750  < — — < 1000) 


Plot  of  Digitized  Velocity  and  Wall  Pressure  Fluctuations  (4950  < — ^ < 5200) 


Digitized  Velocity  and  Wall  Pressure  Fluctuations  (5200  < — - < 5450) 


Fig.  19a  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Velocity  at 

y/6  = .25,  its  VITA  Variance  and  the  Associated  Detector  Function 
tU 
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Fig.  20b  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Velocity  at 
y/6  = its  VITA  Variance,  and  the  Associated  Detector  Function 

(255  < ~ < 505)  68< 


Kig.  20g  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Velocity  at 
y/6  = .14,  its  VITA  Variance,  and  the  Associated  Detector  Function 
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Fig.  22d  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Velocity  at 
y/5  - .025,  its  VITA  Variance,  and  the  Associated  Detector  Function 

(755  < 1005) 
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Fig.  22g  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Velocity  at 
y/6  = .025,  its  VITA  Variance,  and  the  Associated  Detector  Function 


(5205  < ~»<  5455) 
6* 


83< 


Fig.  23b  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Pressure 


at  the  Wall  its  VITA  Variance,  and  the  Associated  Detector  Function 


Fig.  23g  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Pressure 
at  the  Wall,  its  VITA  Variance,  and  the  Associated  Detector  Function 
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23h  Output  of  Digital  Processing  Scheme  Showing  the  Fluctuating  Pressure 
at  the  Wall,  its  VITA  Variance,  and  the  Associated  Detector  Function 
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